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In the dense-neutrino region above the neutrino sphere of a supernova (r < 400 km), neutrino- 
neutrino refraction causes collective flavor transformations. They can lead to "spectral splits" where 
an energy -Egpiit sphts the transformed spectrum sharply into parts of almost pure but different 
flavors. Unless there is an ordinary MSW resonance in the dense-neutrino region, -Bsput is determined 
by flavor-lepton number conservation alone. Spectral splits are created by an adiabatic transition 
between regions of large and small neutrino density. We solve the equations of motion in the 
adiabatic limit explicitly and provide analytic expressions for a generic example. 



PACS numbers: 14.60.Pq, 97.60.Bw 

I. INTRODUCTION 

At large densities, neutrino-neutrino refraction causes 
nonlinear flavor oscillation phenomena with sometimes 
Mrplexing results [l|, 1^, l3|, ij, (Si (fii (tI, [sj, l9|, M, M, M 
ImUdUM- In the region between the neutrino sphere 
and a radius of about 400 km in core-collapse super- 
novae (SNe), the neutrino flavor content evolves dramat- 
ically 0, i, la E [H E El ■ The global features of this 
self-induced transformation are equivalent to the motion 
of a gyroscopic pendulum in flavor space [111, [ll • How- 
ever, this picture does not explain the "spectral splits" 
that have been numerically observed in the transformed 
fluxes [l,[l3,[ll. In a typical case, the primary Ve flux be- 
low a split energy -EgpUt emerges from the dense-neutrino 
region in its original flavor, whereas above i^sput, it is 
completely transformed to (some mixture of and 
Vr), the step at £^spiit being very sharp. (To be specific 
we explore the h'e~Vx system with the atmospheric Am^ 
and the small 13-mixing angle.) 

It has been suggested that an adiabatic transition from 
high to low neutrino density is the primary cause for the 
split 0, [l^ . Dense neutrinos perform synchronized os- 
cillations: all modes oscillate with a common frequency 
Wsynch, even though their individual frequencies vary as 
oj = \Arn? /2E\. Flavor oscillations can be visualized as 
the precession of polarization vectors P^^ in a "flavor B 
field." The "stick together" by the z^-z^-interaction, 
thus forming a collective object that precesses around B. 
The collectivity is lost when the neutrino density de- 
creases. However, if the decrease is slow, all P^^ align 
themselves with or against B in the process of decoupling 
from each other. Eventually they all precess with their 
individual u around B, but without visible consequences 
because of their (anti-)alignment with B. 

We extend this interpretation of the split phenomenon 
in several ways. We (i) show that flavor-lepton number 
conservation determines i?spiit , (ii) solve the equations of 
motion explicitly in the adiabatic limit, and (iii) provide 
an analytic result for a generic case. 



II. EQUATIONS OF MOTION 

We represent the flavor content of an isotropic v-u gas 
by flavor polarization vectors P;j and P;j, where over- 
barred quantities correspond to v. We define their global 
counterparts as P = Jg°° dto and P — duj and 
introduce D = P — P, representing the net lep ton num- 
ber. The equations of motion (EDMs) are [ill, [l^ 

dtP^ = {ujB + AL + /iD) X P„ (1) 

and the same for P(^ with uj ~uj. Here A = V2Gi?ne 
represents the usual matter potential and /i = v^G^nj^ 
the v-v interaction strength, where Ue and n^, are the 
electron and neutrino densities. We work in the mass 
basis where B = (0,0,-1) corresponds to the normal 
and B = (0,0,-fl) to the inverted mass hierarchies. 
The interaction direction L is a unit vector such that 
B • L = cos 26 with 9 being the vacuum mixing angle. 
Unless there is an MSW resonance in the dense-neutrino 
region, one can eliminate AL from Eq. ([T]) by going into 
a rotating frame, at the expense of a small effective mix- 
ing angle [H, [ll|. The only difference for antineutrinos 
is that in vacuum they oscillate "the other way round." 
Therefore, instead of using P^^ we may extend P^^ to 
negative frequencies such that = P > 0) and 

use only P„ with ~oo < uj < +oo. In these terms, 
D — duj Puj, where s^j = sign(cj) = 

After elimination of AL, the EOM for D can be ob- 
tained by integrating Eq. ([T]) with s^: 

/ + 00 
duj s^LuP^ . (2) 
-oo 

It shows that 9t(D • B) = so that = B • D is con- 
served [11]. The in- medium mixing angle above a SN 
core is small and therefore the mass and interaction ba- 
sis almost coincide. Collective effects then only induce 
pair transformations of the form Vei^e i^x^x, whereas 
the excess flux from deleptonization is conserved. 
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III. ADIABATIC SOLUTION 

We rewrite the EOMs in terms of an "effective Hamil- 
tonian" for the individual modes as 

dtF^ = X where B.^ = ujB + /iD. (3) 

In the adiabatic hmit each moves slowly compared to 
the precession of Pj^ so that the latter follows the former. 
We assume that initially all P^^ represent the same flavor 
and thus are aligned. If initially fi is large, every P(^ is 
practically aligned with H^^. Therefore, in the adiabatic 
limit it stays aligned with H^^ for the entire evolution: 



(4) 



which solves the EOMs. Here = \P^\ and H^^ = 
Htj/|Hi^| is a unit vector. Here and henceforth we assume 
an excess flux of neutrinos over antineutrinos, implying 
that initially P^^ and D are collinear and > 0. 

According to Eq. ^ all lie in the plane spanned 
by B and D which we call the "co-rotating plane." In 
the adiabatic limit all P;^, and consequently M, also stay 
in that plane. Therefore we can decompose 



M = bB + ujc'D 
and rewrite the EOM of Eq. ^ as 
dtT) = B X D. 



(5) 



(6) 



Therefore D and the co-rotating plane precess around B 
with the common or "co-rotation frequency" Uc- 

We conclude that the system evolves simultaneously 
in two ways: a fast precession around B determined by 
ujc = Wc{fJ.) and a drift in the co-rotating plane caused 
by the explicit variation. To isolate the latter from 
the former, we go (following Ref. P) into the co-rotating 
frame where the individual Hamiltonians become 



H„ 



{lO — LUc) B + flD . 



(7) 



We use the same notation because the relevant compo- 
nents Hi^z, Hlj±, Dz, and Z)j_ remain invariant. 

Initially (/i — > oo) the oscillations are synchronized, 
— tiJsynch, and all P^^ form a collective P. As ^ de- 
creases, the Pij zenith angles spread out while remaining 
in a single co-rotating plane. In the end — > 0) the 
co-rotation frequency is uj'^ and Eqs. ^ and ([7]) imply 
that all final H^^ and therefore all P^^ with lo > lo^ are 
aligned with B, the others anti-aligned: a spectral split 
is inevitable with cjgpiit = being the split frequency. 
The lengths P^ — \Puj\ are conserved and eventually all 
Ptj point in the ±B directions. Therefore the conserva- 
tion of flavor- lepton number gives us ojspiit, for > 0, 
by virtue of 



P.J duj 



i^split r+oo 

P^duj+ I P^dio. (8) 



In general, ujspht = u;° ^ = Wsynch- 

For individual modes the EOMs given by H^^ are com- 
pletely solved if we find ujc{fi) and D±{fj,), the com- 
ponent transverse to B, since is conserved and 
given by the initial condition. From Eq. Q we infer 
P^±/Pu, = H^±/H^, from Eq. jT]) = ^iD± and 

Hi^z ^ — iUc + ^jlDz so that 



p — 



Integration of the second equation over Si^du gives us 

Pu. 



(9) 
(10) 



r+oo 

1=1 duj Sui — I 
^[{uj-uj^)/^i + D,Y+D\ 



(11) 



Projecting Eq. ([5]) on the ai-y-plane we find lj^ = 
M±/ D±{^) or explicitly 

f^°° dw Suj LO P^x^ duj LO P^±^ 
LOc = — -T-rr = . (12) 



dio s^P^ 



For large /i when the oscillations are synchronized, this 
agrees with the usual expression for Wsynch @5 but it 
changes when the Fi^ spread out in the zenith direction. 
Inserting Eq. ^ into Eq. ^ wc find 



-t-oo 



duj ; 



LoPu, 



V[(I7^7j:)7;it:d:F+^ 



(13) 



Given Dz and a spectrum P^^, we can determine lOc and 
D± from Eqs. (Ill|) and (|13|) for any fi. These equations 
solve the EOMs explicitly in the adiabatic limit. 

We have assumed that all F^ are initially aligned. One 
can relax this restriction and allow some P^^ to have op- 
posite orientation. If different species are emitted from 
a SN core with equal luminosities but different average 
energies, the spectra will cross over so that some range 
of modes is prepared, say, as and another as Vx- 



IV. NEUTRINOS ONLY 

We illustrate the power of our new results with a 
generic neutrino-only example (D = P). The spectrum 
is taken box like with P^ = {2uiq)~^ for < w < 2u}q 
and otherwise. With P^ being conserved we find from 
Eq. ® 



UJQ X 



1 

il-Pz) 



for /i 
for ji 



oo, 
0. 



(14) 



The case = is special because coc = loq remains fixed. 
For Pj = 1 we have wj? = and no flavor evolution. 
We use Pz = 0.5 to show the initial and final Pz{lo) in 
Fig. [T] (left) . The dotted line denotes the adiabatic final 
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FIG. 1: Spectra of polarization vectors (z-compoiient). Thin: initial. Thick: final. Dotted: fully adiabatic. Solid: numerical 
solution as described in the text. Left: Box- like initial v spectrum, large misalignment between B and P, and no i>. Right: 
Box- like v and P spectra, 30% fewer u, small initial misalignment (sin 20 = 0.05), and inverted hierarchy. 
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FIG. 2: Pz,[^(/i) for 51 modes. Left: Box-like ;/-oidy spectrum. Numerical solution of EOMs (top). Analytic adiabatic solution 
{bottom). Right: Box-like i> and u spectra. Numerical solution for i> (top) and u {bottom), here only 6 modes. 
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state where Wgpiit = 0.5 wq- The sohd Hne is from a nu- 
merical solution of the EOMs with ^{t) — fiQ exp(— i/r) 
and T~-^ = 0.03 wq, typical for a SN. We have checked 
numerically that the split indeed becomes sharper with 
increasing t and thus increasing adiabaticity. 

In Fig. [3 we show P^i.zifJ.) for 51 individual modes. 
They start with the common value P^^z = 0.5(2a;o)~^. 
Later they spread and eventually split, some of them ap- 
proaching and the others —1. Some modes first move 
down and then turn around as lOc changes. A few modes 
do not reach ±1 because of imperfect adiabaticity. 

For the box spectrum the integrals Eqs. pT|) and ([Te 
are easily performed and one can extract 



UJc = UJQ + UJqPz 



1 



K e" 
2k 



(15) 



where k = ujq/ij-. For /i ^ cxd and /i — > the limits of 
UJc agree with Eq. (|14p from lepton-number conservation. 
For ^ —> oo we obtain P± = y^l — P^, representing the 
initial condition P = 1, and for /i — s- we find Pi — 0. 

With Eq. ([9]) these results provide analytic solutions 
for the adiabatic Pi^,z(/i). We show examples in Fig. [2] 
(bottom left) for comparison with the numerical solu- 
tion of the EOMs. The agreement is striking and con- 
firms the picture of adiabatic evolution in the co-rotating 
plane. The agreement is poor for modes close to the split 
{lu « w^) at low neutrino densities (/i < ujq) where the 
evolution becomes nonadiabatic. 



V. ADIABATICITY CONDITION 

The speed for the Huj evolution in the co-rotating 
plane is dO^/dt, where cos^^^ = H^±/H^, while P^^ pre- 
cesses with speed H^^. The evolution is adiabatic if the 
adiabaticity parameter = \d9^/dt\H[J^ <^ 1. With 
Eqs. and ^ we find 



UJ - 



D, 



dDi_ 
d^j. 



dujc 
dfi 



3/2 



(16) 



where t„ 



|<iln/i/dt|~i. 
For our neutrino-only 
Eqs. (dni) give dP±/dfi = - 



Pz[l - AK^/ie" - e-'^y] 



P±(Wc 

For 



Pi 



-'± ) box spectrum 
-ujo)/fi^ and dwc/d^ — 
/i ^ Wo we obtain 
dP±/diJL ~ oJq/ij-^ and duc/dfi ^ ^o//^^ ^° thai the last 
term in the numerator of Eq. (|16p dominates: j^j ^ 
P±{uj — ujc)/{hf^^'^). With fi decreasing, 7^ increases and 
at ~ Wo when 7^^ ~ 1, adiabaticity violation begins. 
For n < ujQ the denominator of Eq. (fTB|l gives the depen- 
dence 7^ cx (w — u!c)~^ , and therefore the closer uj to Uc 
the stronger the adiabaticity violation. 



VI. INCLUDING ANTINEUTRINOS. 

As a second generic case we now add antineutrinos. 
One important difference is that even a very small initial 
misalignment between D and B is enough to cause a 
strong effect. Consider a single energy mode for v with 
P = 1 and one for P with P = a < 1 that are initially 
aligned in the flavor direction, now taken very close to 
the mass direction, and assume an inverted hierarchy. 
From the dynamics of the flavor pendulum [ ill , [l^ we 
know that in the end P is antialigned with B, whereas P 
retains a large transverse component because Pz — Pz is 
conserved: The system prepares itself for a spectral split. 

Assuming box spectra for both v and D, we show the 
initial and final Pz,uj in Fig. [1] (right), for the inverted 
hierarchy, sin 20 = 0.05, and a = 0.7. From Eq. one 
infers uj^ — Wgynch = wo (1 + a)/(l — a). For ot = this 
is u)^ — ^ojQ > 2ujo- Therefore, all modes have neg- 
ative frequencies in the co-rotating frame and tilt away 
from B (see also the numerical Puj^z in Fig. [5]). The fi- 
nal split frequency is found from flavor lepton number 
conservation to be Wsput = ^^o (1 ~ Dz + a) ~ ujq 2a, us- 
ing Dz ^ I — a for sin26' <C 1. With a = ^ we find 
^spiit = Y§ in agreement with Fig. [TJ For < a < 1 
we have < Wgput < 2u;o so that the final split always 
occurs among the neutrinos. According to Fig. [2] the split 
starts when the vector D develops a significant transverse 
component, and it proceeds efficiently in a region ~ wo- 

The "wiggles" in the curves in the right panels of F ig.p 
stem from the nutation of the flavor pendulum [ill . [l3 |7 
We have chosen a relatively fast fi{t) evolution (t~^ = 
0.1 Wo), implying poor adiabaticity, to avoid too many 
nutation periods on the plot. For a very slow fj,(t) the 
nutations disappear and the co-rotating frame removes 
the full global evolution of the system. 



VII. DISCUSSION 

We have studied the phenomenon of spectral splits that 
is caused by neutrino-neutrino refraction in the SN dense- 
neutrino region. We have carried previous explanations 
of this novel effect [1, [l2| to the point of explicit solutions 
in the adiabatic limit. 

A spectral split occurs when a neutrino ensemble is 
prepared such that the common direction of the flavor 
polarization vectors deviates from the mass direction. An 
adiabatic density decrease turns all modes below a split 
energy Esput = ATO^/2wspiit into the mass direction, and 
the others in the opposite direction. Remarkably, during 
this phase all modes remain in a single rotating plane, 
even after losing full synchronization. Pspiit is deter- 
mined by lepton number conservation in the mass basis. 

The spectral split is a generic feature of the adiabatic 
evolution when the density changes from large to small 
values. It can appear even in the absense of neutrino- 
neutrino interactions. Indeed, in the usual MSW case 
the evolution to zero density transforms to 1^2 and 
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to v\ for all energies. This corresponds to Wspiit = 0. 
The neutrino-neutrino interactions shift ojspiit to non-zero 
values. 

A spectral split is caused in the SN neutrino (but 
not antineutrino) flux by neutrino-neutrino interactions 
alone, especially during the accretion phase when or- 
dinary MSW resonances occur far outside the dense- 
neutrino region. Later the matter profile may become 
so shallow that the H-resonance moves into this re- 
gion [l^, The simultaneous action of collective 
effects and an ordinary MSW resonance may then cause 
spectral splits for both neutrinos and antineutrinos, lead- 
ing to a rich phenomenology, perhaps modifying r-process 
nucleosynthesis 0, Of course, the fluxes will be fur- 
ther processed by ordinary conversion in the SN enve- 
lope t'^u^ modifying observable signatures. Still, 
observing spectral splits would provide a smoking gun 
signature both for the relevant neutrino properties and, 
if it occurs among antineutrinos at late times, for the oc- 
currence of a shallow density profile above the neutrino 
sphere. 

The neutrino flux emitted by a SN is anisotropic so 
that neutrinos on different trajectories experience dif- 
ferent neutrino-neutrino interaction histories [T3 | that 
would be expected to cause kinematical flavor decoher- 
ence of different angular modes A numerical ex- 



ploration reveals, however, that in a typical SN scenario 
the deleptonization flux suppresses decoherence and the 
evolution is almost identical to that of an isotropic en- 
semble [l^. Our treatment of the spectral evolution is 
apparently applicable in a realistic SN context. 

Collective neutrino oscillation phenomena in a SN may 
well be important for the explosion mechanism, r-process 
nucleosynthesis and may provide detectable signatures in 
a high-statistics signal from the next galactic SN. Build- 
ing on previous ideas, our formalism gives a simple, ele- 
gant and quantitative explanation of seemingly impene- 
trable numerical results. Our approach provides the ba- 
sis for developing a quantitative understanding of realis- 
tic consequences of collective neutrino oscillations for SN 
physics and observational signatures. 



Acknowledgments 

We acknowledge support by the Deutsche Forschungs- 
gemeinschaft (TR 27 "Neutrinos and beyond"), the Eu- 
ropean Union (ILIAS project, contract RII3-CT-2004- 
506222), the Alexander von Humboldt Foundation, and 
The Cluster of Excellence "Origin and Structure of the 
Universe" (Munich and Garching). 



[1] J. Pantaleone, "Neutrino oscillations at high densities," 
Phys. Lett. B 287, 128 (1992). 

[2] S. Samuel, "Neutrino oscillations in dense neutrino gas- 
es," Phys. Rev. D 48, 1462 (1993). 

[3] S. Samuel, "Bimodal coherence in dense selfinteract- 
ing neutrino gases," Phys. Rev. D 53, 5382 (1996) 
hep-ph/ 9604341 . 

[4] Y. Z. Qian and G. M. Fuller, "Matter enhanced anti- 
neutrino flavor transformation and supernova nucleosyn- 
thesis," Phys. Rev. D 52, 656 (1995) astro-ph/9502080 . 

[5] G. M. Fuller and Y. Z. Qian, "Simultaneous flavor trans- 
formation of neutrinos and antineutrinos with dominant 
potentials from neutrino neutrino forward scattering," 
Phys. Rev. D 73, 023004 (2006) astro-ph/0505240 . 

[6] S. Pastor, G. G. Raffelt and D. V. Semikoz, 
"Physics of synchronized neutrino oscillations caused 
by self-interactions," Phys. Rev. D 65, 053011 (2002) 
[hep-ph/0109035 . 

[7] S. Pastor and G. Raffelt, "Flavor oscillations in the su- 
pernova hot bubble region: Nonlinear effects of neu- 
trino background," Phys. Rev. Lett. 89, 191101 (2002) 
[astro- ph/0207281 . 

[8] H. Duan, G. M. Fuller and Y. Z. Qian, "Collective neu- 
trino flavor transformation in supernovae," Phys. Rev. D 
74, 123004 (2006) astro-ph/0511275 . 

[9] H. Duan, G. M. Fuller, J. Carlson and Y. Z. Qian, 
"Simulation of coherent non-linear neutrino flavor trans- 
formation in the supernova environment. L Correlated 
neutrino trajectories," Phys. Rev. D 74, 105014 (2006) 
astro-ph/0606616 . 
[10] H. Duan, G. M. Fuller, J. Carlson and Y. Z. Qian, 



"Coherent development of neutrino flavor in the super- 
nova environment," Phys. Rev. Lett. 97, 241101 (2006) 
astro- ph/0608050 . 

[11] S. Hannestad, G. G. Raffeh, G. Sigl and Y. Y. Y. Wong, 
"Self-induced conversion in dense neutrino gases: Pendu- 
lum in flavour space," Phys. Rev. D 74, 105010 (2006) 
astro- ph/0608695. 

[12] H. Duan, G. M. Fuller, J. Carlson and Y. Z. Qian, "Anah 
ysis of collective neutrino flavor transformation in super- 
novae," jarXiv:astro-ph/0703776 

[13] G. L. FogU, E. Lisi, A. Marrone and A. Mirizzi, "Col- 
lective neutrino flavor transitions in supernovae and the 
role of trajectory averaging," arXiv:0707.1998 [hep-ph]. 

[14] G. G. Raffelt and G. Sigl, "Self-induced decoherence in 
dense neutrino gases," Phys. Rev. D 75, 083002 (2007) 
hep-ph/0701182 . 

[15] A. Esteban-Pretel, S. Pastor, R. Tomas, G. G. Raf- 
felt and G. Sigl, "Decoherence in supernova neu- 
trino transformations suppressed by deleptonization," 
arXiv:0706.2498 [astro-ph]. 

[16] G. Sigl and G. Raffelt, "General kinetic description of 
relativistic mixed neutrinos," Nucl. Phys. B 406, 423 
(1993). 

[17] A. S. Dighe and A. Y. Smirnov, "Identifying the neutrino 
mass spectrum from the neutrino burst from a super- 
nova," Phys. Rev. D 62, 033007 (2000) hep-ph/9907423; . 

[18] A. Dighe, "Supernova neutrinos: Production, propaga- 
tion and oscillations," Nucl. Phys. Proc. Suppl. 143, 449 
(2005) hep-ph/0409268| . 



